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Abstract 



The breaking of supersymmetry due to singular potentials in supersymmetric quantum me- 
chanics is critically analyzed. It is shown that, when properly regularized, these potentials 
respect supersymmetry, even when the regularization parameter is removed. 
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1 Introduction 



Supersymmetry is a beautiful and, simultaneously, a tantalizing symmetry [1-7]. On the one hand, 
supersymmetry leads to field theories and string theories with exceptional properties [8-9]. The im- 
proved ultraviolet behavior, the natural solution of the hierarchy problem are just a few of the nice 
features of supersymmetric theories. On the other hand, supersymmetry also predicts degenerate 
superpartner states (superpartner states with degenerate mass) corresponding to every physical 
particle state of the theory. Unfortunately, such superpartners with degenerate masses are not ob- 
served experimentally and, consequently, one expects that supersymmetry must be spontaneously 
(dynamically) broken much like the spontaneous breaking of ordinary symmetries in physical the- 
ories. However, unlike ordinary symmetries, spontaneous breaking of supersymmetry has so far 
proved extremely difficult in the conventional framework. Consequently, in the context of super- 
symmetry, one constantly looks for alternate, unconventional methods of breaking of this symmetry 
[6-7]. There is, of course, the breaking of supersymmetry due to instanton effects which is well un- 
derstood. However, several authors, in recent years have suggested that supersymmetry may be 
broken in the presence of singular potentials or boundaries in a nonstandard manner [10-12]. We 
will explain the details of the proposed mechanism later, but the gist of the argument is that in such 
systems, the superpartner states may not belong to the physical Hilbert space thereby leading to a 
breaking of supersymmetry. This would, of course, explain why the superpartner states would not 
be observable. Even more interesting is the possibility that since such a breaking is nonstandard, 
the usual theorems of supersymmetry breaking may not apply and the ground state may continue 
to have zero energy thereby leading to a solution of the cosmological constant problem also. 

The examples, where such a breaking has been discussed, are simple quantum mechanical models 
which nonetheless arise from the non-relativistic limit of some field theories. It is for this reason 
that, in an earlier paper, we had examined [13] a candidate relativistic 2 + 1 dimensional field 
theory to see if the manifestation of such a mechanism was possible in a field theory. However, 
a careful examination of the theory revealed that supersymmetry prevails at the end although it 
might appear naively, in the beginning, that supersymmetry would be broken in the nonstandard 
manner. This prompted us to re-analyze the quantum mechanical models, where this mechanism 
was demonstrated, more carefully in order to understand if supersymmetry is truly broken in 
these models and if so what may be the distinguishing features in these models from a relativistic 
field theory. A systematic and critical examination, once again, reveals that when carefully done, 
supersymmetry is manifest even in such singular quantum mechanical models which is the main 
result of this paper. 

Since our discussion would be entirely within the context of one dimensional supersymmetric 
quantum mechanics, let us establish the essential notations here. Given a superpotential, W(x), 
we can define a pair of supersymmetric potentials as 



where "prime" denotes differentiation with respect to x. With h = 1 and m = 1, we can, then, 




(1) 
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define a pair of Hamiltonians which describe a super symmetric system as 

1 d 2 

1 d 2 



H- = o+V- ( 2 ) 

2dx 2 K ' 



In fact, defining the supercharges as 



1 U.wmV ot = -L( d 



we recognize that we can write the pair of Hamiltonians in eq. (§) also as 

H+ = H- = QQ ] (4) 



It is clear now that if l^— ) represents an eigenstate of H— with a nonzero energy E, then, 
Qt | ■)/>_) would be an eigenstate of H + with the same energy, namely, 

= QQ^-)=E\^) 
F + (Qt|^_)) = Qt Q (Qt| V ,_)) = Q\{QQ^_)) = E(Q^)) (5) 

In other words, the states \ip-} and Q^\ip-} (alternately, we can also denote them as QI^M-) and 
\ip+) respectively) would correspond to the degenerate superpartner states. All the eigenstates of 
the two Hamiltonians H + and H- would be degenerate except for the ground state with vanishing 
energy which would correspond to the state satisfying 

QhM=0, or, g f |V-) = (6) 



For a given superpotential, at most one of the two conditions in eq. (|6|) can be satisfied (that 
is, at most, only one of the two conditions in (|6|) would give a normalizable state). Namely, the 
ground state with vanishing energy is unpaired and can belong to the spectrum of either H + or H_ 
depending on which of the conditions leads to a normalizable state. This corresponds to the case 
of unbroken supersymmetry. For example, with W(a;) = —lux, the pair of Hamiltonians 

describe the supersymmetric harmonic oscillator and it can be easily checked that the normalizable 
ground state belongs to the spectrum of H + (In fact, we will choose H + to have the ground state 
throughout this paper.). 

If, on the other hand, the superpotential is such that neither of the states in eq. (||) is normal- 
izable, then, supersymmetry is known to be broken by instanton effects [6]. For example, this can 
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happen when the superpotential is an even polynomial. In contrast, the new mechanism described, 
in the presence of a boundary or a singular potential [10-12], corresponds to the case where the 
action of the supercharges takes a state out of the physical Hilbert space (although QQ^ and Q^Q 
belong to the Hilbert space) so that, say, if |Y>+) is in the Hilbert space, Q\ip+) would not belong to 
the Hilbert space leading to the breaking of the degeneracy of states and, therefore, supersymmetry. 

In this paper, we carefully analyze the models [10-12] where this mechanism is thought to 
be operative. We note that a quantum mechanical potential with a singular structure is best 
studied with a regularization because this brings out the correct boundary conditions naturally. 
Furthermore, the regularization must be chosen carefully preserving supersymmetry when one is 
dealing with a supersymmetric potential with a singular structure. Keeping this in mind, we 
examine a supersymmetric model with a boundary, namely, the supersymmetric oscillator on the 
half line [10] in section 2. We briefly recapitulate the results of a harmonic oscillator on the half line 
and then show through a careful analysis that supersymmetry is, in fact, manifest in such a model 
in spite of the boundary. In section 3, we analyze a model with a singular potential, namely, the 
oscillator with a 4j potential. This potential is quite interesting and, in spite of several studies in 
the literature [14-16], the analysis is not quite complete and we present a systematic and complete 
analysis of this problem. In section 4, we take up the study of the supersymmetric oscillator with a 
4j potential [11] and show that supersymmetry is, in fact, manifest in this system as well, in spite 
of the singular nature of the potential. There are several interesting aspects of this model which 
emerge from a careful analysis which we bring out. In section 5, we solve this problem algebraically 
as well which supports the analysis of section 4. In section 6, we describe the solution to a puzzle 
raised in the literature [11] in the context of the supersymmetric oscillator with a 4j potential and 
present our conclusions in section 7. 



2 Super "Half" Oscillator 



To understand the super "half" oscillator, it is useful to recapitulate briefly the results of the "half" 
oscillator. Let us consider a particle moving in the potential 

V(x) = l U^ x2 -") for *>° (8) 
w 1 oo for x < v ' 

The spectrum of this potential is quite clear intuitively. Namely, because of the infinite barrier, we 
expect the wave function to vanish at the origin leading to the conclusion that, of all the solutions 
of the oscillator on the full line, only the odd solutions (of course, on the "half" line there is no 
notion of even and odd) would survive in this case. While this is quite obvious, let us analyze the 
problem systematically for later purpose. 

First, let us note that singular potentials are best studied in a regularized manner because this 
is the only way that appropriate boundary conditions can be determined correctly. Therefore, let 
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us consider the particle moving in the regularized potential 

\{uj 2 x 2 — u) for x > 



V(x) = (9) 



\ for x < 

with the understanding that the limit |c| — > oo is to be taken at the end. The Schrodinger equation 

(-5^ + ^)) ^)=*) 

can now be solved in the two regions. Since |c| — * oo at the end, for any finite energy solution, we 
have the asymptotically damped solution, for x < 0, 

^ n \x) = Ae^- 2 ^ x (10) 



Since the system no longer has reflection symmetry, the solutions, in the region x > 0, cannot 
be classified into even and odd solutions. Rather, the normalizable (physical) solution would 
correspond to one which vanishes asymptotically. The solutions of the Schrodinger equation, in the 
region x > 0, are known as the parabolic cylinder functions [17] and the asymptotically damped 
physical solution is given by 

^\x) =BU(-(- + h,V2^x) (11) 

The parabolic cylinder function, U(a, x), of course, vanishes for large values of x. For small values 
of x, it satisfies 



U(a,x) x -° 



2 I(2a+l) r( 3 + «> 



tf'(a,x) ^ -— r p- (12) 



It is now straightforward to match the solutions in eqs. ( |10j , |11| ) and their first derivatives across 
the boundary at x = and their ratio gives 



i _ i r(-£) 



It is clear, then, that as |c| — > oo, this can be satisfied only if 



(13) 



-i + 5 ICj ^-"' " = 0,1,2,... (14) 

In other words, when the regularization is removed, the energy levels that survive are the odd ones, 
namely, (remember that the zero point energy is already subtracted out in (||) or (^)) 

e n = w(2n + l) (15) 
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The corresponding physical wave functions are nontrivial only on the half line x > and have the 
form 

1> n {x) = B n U(-(2n + |), v^x) = B n e"^ 2 H 2n+1 (y/^x) (16) 

Namely, only the odd Hermite polynomials survive leading to the fact that the wave function van- 
ishes at x = 0. Thus, we see that the correct boundary condition naturally arises from regularizing 
the singular potential and studying the problem systematically. 



We now turn to the analysis of the supersymmetric oscillator on the half line. One can define 
a superpotential [10] 

W(x) = { ~ UX 5° r X > ° (17) 

v ' I oo for x < v ' 

which would, naively, lead to the pair of potentials 

| iCuVTuO for x>0 

1 oo for x < v ' 

Since, this involves singular potentials, we can study it, as before, by regularizing the singular 
potentials as 



V+(x) 



\(u) 2 x 2 — oS) for x > 



-y- for x < 

^(oj 2 x 2 + w) for x > 



2 

with the understanding that \c±\ — > oo at the end 



V-(x) = { n (19) 

for x < 



The earlier analysis can now be repeated for the pair of potentials in eq. (|l9|) . It is straightfor- 
ward and without going into details, let us simply note the results, namely, that, in this case, we 
obtain 

e+,n = uj(2n + l) tp+, n (x) = B+^e"^ H 2n +i(Vux) 

e_, n = 2w(n + l) tp-,n(x) = B^ n e~-^ H 2n+1 {^x) (20) 

Here n = 0, 1, 2, • • •. There are several things to note from this analysis. First, only the odd Hermite 
polynomials survive as physical solutions since the wave function has to vanish at the origin. This 
boundary condition arises from a systematic study involving a regularized potential. Second, the 
energy levels for the supersymmetric pair of Hamiltonians are no longer degenerate. Furthermore, 
the state with e = no longer belongs to the Hilbert space (since it corresponds to an even Hermite 
polynomial solution). This leads to the conventional conclusion that supersymmetry is broken in 
such a case and let us note, in particular, that in such a case, it would appear that the superpartner 
states do not belong to the physical Hilbert space (Namely, in this case, the supercharge is an odd 
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operator and hence connects even and odd Hermite polynomials. However, the boundary condition 
selects out only odd Hermite polynomials as belonging to the physical Hilbert space.). 

There is absolutely no doubt that supersymmetry is broken in this case. The question that 
needs to be addressed is whether it is a dynamical property of the system or an artifact of the 
regularization (and, hence the boundary condition) used. The answer is quite obvious, namely, 
that supersymmetry is broken mainly because the regularization (and, therefore, the boundary 
condition) breaks supersymmetry. In other words, for any value of the regularizing parameters, c± 
(even if |c+| = |c_|), the pair of potentials in eq. (|l^) do not define a supersymmetric system and 
hence the regularization itself breaks supersymmetry. Consequently, the breaking of supersymmetry 
that results when the regularization is removed cannot be trusted as a dynamical effect. 



Regularized Superpotential 



Another way to understand this is to note that for a supersymmetric system, it is not the potential 
that is fundamental. Rather, it is the superpotential which gives the pair of supersymmetric poten- 
tials through Riccati type relations. It is natural, therefore, to regularize the superpotential which 
would automatically lead to a pair of regularized potentials which would be supersymmetric for any 
value of the regularization parameter. Namely, such a regularization will respect supersymmetry 
and, with such a regularization, it is, then, meaningful to ask if supersymmetry is broken when the 
regularization parameter is removed at the end. With this in mind, let us look at the regularized 
superpotential 

W(x) = -ujx6(x) + cd(-x) (21) 

Here c is the regularization parameter and we are supposed to take \c\ — > oo at the end. Note that 
although, at this level, both the signs of the regularization parameter are allowed, existence of a 
normalizable ground state (see eq. @) selects out c > (otherwise, the regularization would have 
broken supersymmetry through instanton effects as we have mentioned earlier) . 



The regularized superpotential now leads to the pair of regularized supersymmetric potentials 

1 r. 



V-(x) 



{to 2 x 2 - uj)6(x) + c 2 6(-x) - c5(x) 
{u?x 2 + oj)0(x) + c 2 6(-x) + c6(x) 



(22) 



which are supersymmetric for any c > 0. Let us note that the difference here from the earlier case 
where the potentials were directly regularized (see eq. fll9|)) lies only in the presence of the 8(x) 
terms in the potentials. Consequently, the earlier solutions in the regions x > and x < continue 
to hold. However, the matching conditions are now different because of the delta function terms. 
Carefully matching the wave function and the discontinuity of the first derivative across x = for 
each of the wavefunctions and taking their ratio, we obtain the two conditions 

1 1 r(- 



2u)> 



(c 2 -2e+) 1 /2 



2 v^r(-|± + i 



(23) 
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(24) 




e +,n 



= 2um 



n = 0, 1,2, ••• 



,n 



2w(n + 1) 



(25) 



The corresponding wave functions, in this case, have the forms 



(26) 



This is indeed quite interesting for it shows that the spectrum of H + contains the ground state 
with vanishing energy. Furthermore, all the other states of H + and H- are degenerate in energy 
corresponding to even and odd Hermite polynomials as one would expect from superpartner states. 
Consequently, it is quite clear that if the supersymmetric "half" oscillator is defined carefully 
by regularizing the superpotential, then, supersymmetry is manifest in the limit of removing the 
regularization. This should be contrasted with the general belief that supersymmetry is broken in 
this system (which is a consequence of using boundary conditions or, equivalently, of regularizing 
the potentials in a manner which violates supersymmetry). 



Alternate Regularization 

Of course, we should worry at this point as to how regularization independent our conclusion really 
is. Namely, our results appear to follow from the matching conditions in the presence of singular 
delta potential terms and, consequently, it is worth investigating whether our conclusions would 
continue to hold with an alternate regularization of the superpotential which would not introduce 
such singular terms to the potentials. With this in mind, let us choose a regularized superpotential 
of the form 



Here A is the regularization parameter and we are to take the limit |A| — > oo at the end. Once 
again, we note that, although both signs of A appear to be allowed, existence of a normalizable 



This regularized superpotential would now lead to the pair of supersymmetric potentials of the 
form 




(27) 





(28) 
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There are no singular delta potential terms with this regularization. In fact, the regularization 
merely introduces a supersymmetric pair of oscillators for x < whose frequency is to be taken to 
infinity at the end. 



Since there is a harmonic oscillator potential for both x > and x < 0, the solutions are 
straightforward. They are the parabolic cylinder functions which we have mentioned earlier. Now 
matching the wave function and its first derivative at x = for each of the Hamiltonians and taking 
the ratio, we obtain 



VXr(-^ + i) v£r(-£ + i) 

It is clear now that, as A — > oo, eqs. ( |29| ) and (pQ) give respectively 

e+ )n = 2um 



The corresponding wave functions are given by 



1>-, n {x) = B_ :n e-^ x2 H 2n+1 (^x) (32) 

These are, of course, the same energy levels and wave functions as obtained in eqs. ( p5| ) and 
(Pq) respectively showing again that supersymmetry is manifest. Furthermore, this shows that 
this conclusion is independent of the regularization used as long as the regularization preserves 
supersymmetry which can be achieved by properly regularizing the superpotential. 



(29) 
(30) 



n = 0, 1,2, ••• 

2w(n + l) (31) 



3 Oscillator with \ Potential 



In the last section, we showed that, in the presence of one kind of singularity, namely, a boundary, 
supersymmetry is unbroken. In what follows, we will study another class of supersymmetric models, 
namely, the supersymmetric oscillator with a -\ potential, where there is a genuine singularity in 
the potential not necessarily arising from a boundary. A naive analysis of this model [11] also 
shows that supersymmetry is broken by such a singular potential (for certain parameter ranges). 
However, this conclusion can be understood, again, as a consequence of regularizing the potential 
which, as we have seen before, does not respect supersymmetry. In stead, we will show through 
a careful analysis that, when the superpotential is regularized, supersymmetry is manifest in this 
model as well (with a lot of interesting features). In this section, however, we will systematically 
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analyze only the quantum mechanical system corresponding to an oscillator in the presence of a \ 
potential (postponing the discussion of the supersymmetric case to the next section). This system 
has been analyzed by several people [14-16] and the most complete analysis appears to be in ref. 
[16]. However, we feel that, while the energy levels derived in [16] are correct, the wave functions 
are not (namely, the extensions of the solutions from the positive to the negative axis are incomplete 
and the wave functions, of course, become quite crucial when one wants to extend the analysis to a 
supersymmetric system) and, consequently, we present a careful analysis of this system regularizing 
the singular potential in a systematic manner. With the supersymmetric system in mind (to follow 
in the next section), we write the potential for the system as (with h = m = uj = 1) 



V{x) 



g(g + i) , 2 

n h X 



2g + l 



Consequently, the Schrodinger equation that we want to study has the form 

" d 2 9(9 + 1) 



dx 2 



x 2 + (2e + 20-1) 



ip(x) 







(33) 



(34) 



The singular potential is repulsive for g > or g < — 1 while it is attractive for — 1 < g < 0. 
Furthermore, for ease of comparison, let us make the identifications with the notations of ref. [16] 
(note that our energies are shifted since we have in mind the supersymmetric system to study later) 

2 V 4 7 2 
E = e + 9-\ (35) 



It is also worth noting here that the Schrodinger equation in (34) is invariant under 



g -(g + 1) 

e <-> e + 2g + l (36) 

This symmetry, of course, would also be reflected in the solutions. Furthermore, the fixed point 
of this symmetry, namely, g = — | separates the two branches (namely, for every value of A there 
exist two distinct values of g corresponding to two distinct branches separated at the branch point) 
in the parameter space. 



Regularized Potential 



The Schrodinger equation in fl34| ) can be solved quite easily for x > as was also done in [16]. 
However, to determine correctly how this wavefunction should be extended to the negative axis, it 
is more suitable to regularize the potential near the origin and study the problem carefully. Let us 
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consider a potential of the form 



V(x) 



"g(g+i) 



+ x 2 -25 + l 



S(9±l) +R 2_ 2g+ l 



for Ixl > R 



for Ixl < i? 



(37) 



Namely, we have regularized the potential in a continuous manner preserving the symmetry in 
eq. ( |36| ) with the understanding that the regularization parameter R —* at the end. With this 
regularization, the Schrodinger equation has to be analyzed in three distinct regions. However, 
since the potential has reflection symmetry, we need to analyze the solutions only in the regions 
—R < x < R and x > R. 



The potential is a constant in the region —R < x < R and hence the Schrodinger equation is 
quite simple here. The solutions can be classified into even and odd ones and take the forms 

ijj( H > ven (x) = A(R) cosh kx 

^(ll)odd^ = B (fy s - mhKX ( 38 ) 

where we have defined 

«=^ + fl J -(2e + 2 S -l)«^H (39 ) 

Since R is small (and we are to take the vanishing limit at the end), the last equality holds only 
if g ^ Oor — 1 which we will assume. The special values of g corresponding to the absence of a 
singular potential have to be treated separately and we will come back to this at the end of this 
section. We note here that the normalization constants, A and B, can, in principle depend on the 
regularization parameter which we have allowed for in writing down the form of the solutions in 
eq. ©. 

The potential is much more complicated in the region x > R. However, if we make the definitions 

f = x 2 , V (7) = e'^ ft u(£) (40) 

where 

s = (g + 1) or - g (41) 
the Schrodinger equation of (|4|) takes the form 

Su . 1 .du I . . . , . 

^ + ( S +2-fc + 2 (e + 9 - S - 1)u = ° (42) 



An equation of the form 



d 2 u ., 







du 



au 



(43) 
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is known as the confluent hyper geometric equation [17] and the only solution of this equation which 
is damped for large values of the coordinate has the form 

U ("> 6 ' = vhi~ b \\ M ( a > 6 ' ® + ^fr-t^MO. + a-b,2-b,0 (44) 
r(l + a — b) T{a) 

Here M(a, b, £) are known as the confluent hypergeometric functions which have the series expansion 



and satisfy 



Af(a, 6,0^1 (45) 



It is clear now that eq. (42) simply is the confluent hypergeometric equation and the asymptot- 
ically damped physical solutions are nothing other than the U (a, b, £) functions with appropriate 
parameters. Since s takes two possible values (see eq. (glD), it would seem that there would be 



two independent solutions of eq. (42). However, it can be easily checked that the two solutions 
corresponding to the two values of s are really proportional to each other and not independent. 
Thus, we can write the general solution of the Schrodinger equation, for x > 0, as 

^W( x ) = C(R)e-& (46) 
rb-^c+i,^ e _, 3 ^ , Tig + $ _ e _,1 



T{\-g-^f ^ 2' y ' 2'" ; ' r(l-f)"" MK 2 " 2' y ' 2 

Once again, we have allowed for a dependence of the normalization constant, C, on the regulariza- 
tion parameter, R. However, for a nontrivial solution to exist, we require that 

C(R) ^°C^0 



S+ 1 M(l - * g + 2, x 2 ) + M(i - 5 - - 5 + i, x 2 ) 



So far, we have the general solutions, in the two regions, where energy is not quantized and 
which should arise from the matching conditions. Furthermore, we have not bothered to evaluate 
the solution in the region x < —R which clearly would be the same as in the region x > R. However, 
the matching conditions would determine how we should extend the solutions in the region x > R 
to the region x < —R. Therefore, let us now examine the matching conditions systematically since 
there are two possible cases. 

(i) Even Solution 

We can match the even solution of the region —R<x<R and its derivative with those of the 
region x > R at x = R. Taking the ratio and remembering that R is small (which is to be taken to 
zero at the end), we obtain to the leading order in R 

, (g + V&^BP+i-g^BT' 

g(g + 1) tanh ^ g{g + 1) = % pi 2> (47) 
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Since the left hand side is independent of R, for consistency, the right hand side must also be and 
this can happen in two different ways. 



First, for g > — ^, it is clear that relation (fiTl) can be satisfied if (we assume from now on that 
n = 0,l,2,---.) 

l-~ = -n + h(g)R 2 9 +1 
or, e n = 2(n + l)-2f 1 (g)R 2 a+ 1 (48) 
with a suitable choice of fi(g). 

On the other hand, for g < — i, if 



or, e n = (2n-2g + l)-2f 2 (g)R- 29 - 1 (49) 



relation fl47| ) can be satisfied with a suitable choice of /2(g)- It is clear that the two possible 
branches of the solution simply reflect the symmetry in eq. p6|) . 



This analysis shows that when the regularization is removed (namely, R — > 0), we have an even 
extension of the solution of the forms 

9>-\: e n = 2(n + l) (50) 

with 

Mx)=Cn r(-g-l-n) e 2 M ^a + - 2 ,x 2 )\ lx]g+1 for x<Q (51) 



ff<~2 : e n = 2n-2 5 + l (52) 

with 

, / x ^ ^(5+9) _i T 2 1 o N I for x > 

^ w (z) = C ra , 1 N e 2" M(-n, -0+ -,x ) < , , „ (53) 

v ; T(g + |-n) 2' J \ \x\ 9 for x < v ' 



(ii) Odd Solution 

We can also match the odd solution of the region —R < x < R and its derivative with those of 
the region x > R at x = R and taking the ratio, we obtain to leading order 



g(g + 1) coth ^g(g + 1) = % ^ x \ (54) 
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Clearly, the analysis following from eq. (47) goes through identically so that we conclude that in 
the limit R — ► 0, we have an odd extension of the solution of the forms 



9> ~r 



e n = 2(n + 1) 



(55) 



with 



ifj n (x) = C n 



T{-g-\-n 



■e~2 x M(-n,g + -,x 2 



r-9+l 



for x > 
for x < 



(56) 



9< ~2 : 



with 



i/j n (x) = C n 



T{g + \-n) ( 



1„2 
2 



M(- 



•n, 



2n - 2# + 1 



" 5+ 2' X 



for x > 
for x < 



(57) 
(58) 



Understanding of the Result 



The conclusion following from this analysis, therefore, is that every energy level of this system is 
doubly degenerate. Both even and odd extensions of the solution are possible for every value of the 
energy level. The energy levels, as given in eqs. ( |50|) and ( |52| ) (or, alternately, ( |55| ) and (|57|)) are, 
of course, identical to those obtained in [16]. The crucial difference is in the structure of the wave 
functions, namely, that both even and odd extensions of the solution are possible for every value 
of the energy (Incidentally, the solutions we have obtained in terms of confluent hyper geometric 
functions also coincide with generalized Laguerre polynomials as was obtained in ref. [16].). It is 
crucial, therefore, to ask if such a conclusion is physically plausible. To understand this question, 
let us recapitulate the results from a simple quantum mechanical model which is well studied. 
Namely, let us look at a particle moving in a potential of the form 



V(x) 



■jS(x) for \x\ < a 
oo for Ixl > a 



It is well known that the solutions of this system can be classified into even and odd ones with 
energy levels (h = m = 1) 



n 2 7r 2 



2(a + i) 2 



n 2 7r 2 



2a 2 



The even and the odd solutions, of course, have distinct energy values for any finite strength of the 
delta potential. However, when 7 — > 00, both the even and the odd solutions become degenerate in 
energy. Namely, a delta potential with an infinite strength leads to a double degeneracy of every 
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energy level corresponding to both even and odd solutions. The connection of this example with 
the problem we are studying is intuitively clear. Namely, we can think of 

£k±il = Km sto+i) = lim fl_^_) 

x z t]— >o x A + r\ z r)-^o \ rj J \ir x z + r/ z J 

It is clear that for g Oor — 1, the singular -\ potential behaves like a delta potential with an 
infinite strength and it is quite natural, therefore, that this system has both even and odd solutions 
degenerate in energy. 



It is also clear from this analysis that it is meaningless to take the g = Oor — 1 limit from the 
results obtained so far simply because the characters of the two problems are quite different. As 
we have argued, for any finite value of g not coinciding with those special values, the potential 
behaves, at the origin, like a delta potential of infinite strength while for the special values, there 
is no such potential. The two cases are related in a drastically discontinuous manner. As a result, 
one cannot treat the 9 ^ 9 ^^ as a perturbation and obtain the full, correct solution simply because 
there is nothing perturbative (small) about this potential for any "nontrivial" value of g. Another 
way of saying this is to re-emphasize what we have already observed following eq. (|39|), namely, 
the character of k and, therefore, the matching conditions change depending on whether or not g 
differs from the special values 0, —1. 



To see how the standard results of the harmonic oscillator would emerge from this analysis, let 
us work out the case only for g = 0. In this case, the matching conditions for the even solution 
would lead to the relation for the ratios 

r /i_ 2 i\ ~~ (2 £ ~~ 1)-^ r(i-i) 
-(*-'>« = 4tF|) T rtt, < 59 > 

which can be satisfied only if 

~ - \ = ~n, or, e n = (2n + 1) (60) 

We recognize these to be the even levels of the oscillator (remember the shifted zero point energy 
of the system in eq. (|33|) for g = 0). Similarly, matching the odd solution and its derivative leads 
to 



which can be satisfied only if 



r( i- f ) t rlT^f) 



l-| = -n, or, e n = 2(n + l) (62) 



These are, of course, the odd energy levels of the harmonic oscillator. (The corresponding wavefunc- 
tions are the even and odd Hermite polynomials respectively.) There is no longer any degeneracy 
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of these levels. In other words, the matching conditions change drastically and so do the solutions 
of the problem depending on whether or not g equals one of the special values 0,-1 and the results 
for the special values cannot be (and, in fact, should not be) obtained from the general result in a 
limiting manner. 



4 Supersymmetric Oscillator with \ Potential 



In this section, we will analyze the supersymmetric version of the case studied in the last section. 
Let us consider a superpotential of the form [11] 



W(x) = --x 

X 

so that the pair of supersymmetric potentials would have the form 

V + (x) 
V-{x) 



(63) 



9{9 ~ l) + x*-2g-l 



x* 



+ x-2g + l 



(64) 



To analyze this problem, we should, of course, regularize the superpotential. However, even before 
introducing the regularization, let us observe some general features associated with this system, 
namely, that this supersymmetric system would have a ground state satisfying 



or, ip (x) 







- I e 

a 



(65) 



which would be damped for asymptotically large values of the coordinate. On the other hand, from 
the behavior of this wave function for small values of x, it is clear that normalizability of the wave 
function would require that g > — ^. We also note that since the supersymmetric system involves 
two Hamiltonians with different g dependence, the symmetry observed in the previous section, 
namely, g <-> —(g + 1) (or, alternately, g <-> —(g — 1)) cannot be a symmetry of the whole system 
(Another way of saying this is to note that the superpotential has no such symmetry. We will come 
back to this question later in this section.). Furthermore, we cannot naively take over the results 
from the previous section since, as we have seen earlier, regularizing the potential directly may not 
respect supersymmetry. 



Therefore, to study this problem systematically, we regularize the superpotential as 

W(x) = 9(x - \R\) (f - *) + 0{R - \x\) (| " ^) | (66) 

Here, as before, R is the regularization parameter which should be taken to zero at the end and 
we have regularized the superpotential such that it is continuous across the boundary. This has 
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the nice feature that there are no delta potential terms in the potentials. In fact, the regularized 
superpotential leads to the pair of super symmetric potentials of the forms 



V+(x) = \ 



9{x - \R\) + x 2 - 2g - l) + 9{R - |x|) ((X _ i) x 2 + - ^ 

'C* - \R\) (^1) +x >- 2g + l ) + e iR - |x|) (U - 1) 2 x 2 - (J- - l)) 



(67) 



It is worth noting here that for g = 0, the singular potential at the origin is not present in both 
the Hamiltonians (namely, it is truly not there). On the other hand, for g = lor — 1, although the 
singular potential disappears from only one of the Hamiltonians, the complete system remembers 
about it through supersymmetry as is clear from the structure of the regularized potentials. 



Spectrum of 



Let us now analyze the spectrum of the two different Hamiltonians systematically. First, let us 
note that, for the Hamiltonian H + , the potential is that of a harmonic oscillator in the region 
—R < x < R. We can write down the even and the odd solutions in this region as [17] (We will 
assume throughout that g ^ 0) 



^° d \x) = B + (R)xe~^> 2 M(-f^ + l, 3 -,(j- 2 -l)x^) (68) 



The solution in the region x > R can also be obtained from an analysis as given in the earlier 
section and leads to 

•$){x) = C + (R)e~^ x ' 2 (69) 
T{ ~ 9 + h) -x 9 M(- e -±,g + 1, x 2 ) + r if~^ -^M(l - g - e -±, -g + I x 2 ^ 



Y{\- g -'-±) y 2' y 2' ; r(-^) y 2 y 2 ' 2 ' 

Once again, we can match the solutions in the two regions and their derivatives across x = R to 
obtain the relevant quantization conditions. As before, there are two possibilities. 

(i) Even Solution 

If we match the even solution in eq. ( |68[ ) and its derivative with those of eq. ( |69] ) at x = R and 
take the ratio, we obtain to leading order in R (remember it! is small and is to be taken to zero at 
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the end), 



r(-g+|) R9 _ 1/ _ i) r <J>-?) R - g +i 

= i-l^lil £E1T (70) 

r(*-*-*)^ + WW) R 
Clearly, this can be satisfied only if (n = 0, 1, 2, • • •) 

- y = -n, or, e +in = 2n (71) 

Thus, we obtain that there exists an even extension of the solution of the form (when the regular- 
ization is removed) 

/ \ ^ r(— g + n) _i T 2 1 2 N I for x > ,„> 

^,»(x) = C7 + ,. r( _^ + i^ w) e >■ M(-n, 9 + -,x 2 )| ^ for x<Q (72) 

with the energy levels given by 

e +5 „ = 2n (73) 



Furthermore, we note from equation (|72| ) that normalizability of the solution restricts that this 
solution is physical only for g > — |. 

(ii) Odd Solution 



We can also match the odd solution in eq. ( |6q ) and its derivative with those of eq. (^) and the 
ratio leads to the appropriate quantization conditions. Without going into technical details, let us 
simply note here that there are two possibilities in this case. For g > |, we have in the limit R — > 

e+, n = 2n (74) 

with 

/ \ ^ r(— 5+,) _i j , ,, 1 9 . I for x > 

^) =g ^ r(- g +^-n) e M (-»^+2'^{ -|x| s for x<0 < 75 > 

On the other hand, for g < i, 

e+, n = 2n - 2 5 + 1 (76) 

with 

, ^ r (5-i) -1.3,,, 3 2 , f x~9 +1 for x>0 

^^ = C+ - r (g -l-n) e 2 M (-»'-2+2'*){ -|x|-^ for x<0 ^ 

This completes the analysis of the spectrum for H + . We see that it consists of a set of even 
solutions with vanishing energy as we would expect from supersymmetry, but it also contains 
additional physical solutions. 
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Spectrum of H- 



We can similarly analyze the spectrum for the super symmetric partner Hamiltonian H-. The even 
and the odd solutions in the region —R < x < R have the forms 

(x) = A-(«)^^ l2 ^,i(|-V) 

^ I)odd{x) = B _ {R)xe -^-^ M (-^ + y-,(^-iy) (78) 

while the solution in the region x > R now takes the form 

^\ x ) = C_(i?)e-^ 2 (79) 



The matching of the solutions and their derivatives can now be done in the standard manner 
and we simply give the results in the limit R — > here. There exists a normalizable even solution 
only when g < i of the form 

/ c ^ n r (#+l) -\x 2 , 1 2x / for ^ > , Qn , 

w_ n (x) = G_ n — — -e 2 X M(— n, — g + -, a; ) < . ._„ . _ (80) 

' ' T(p + i-n) 2 [ \x\ 9 for x < v ; 

with 

e_ „ = 2n - 2g + 1 (81) 
On the other hand, there are two possible odd solutions. For g > — |, we have 

/ / \ /-< r ( _ 3 ~ |) _I x a ,,, 3 2 J x s+1 for x>0 , . 

tb- n (x) = C- n — e a 31 M(—n, g + -, ar) < , .„,-, . _ (82) 

Y ,nv ; ' n r(-c/-i-n) 2 ; | -\x\9+ l for x < y ' 

with 

e_, n = 2(n + 1) (83) 

while, for g < — |, we have 

/ / \ r(^+i) _i„2 1 2 | x~ 9 for x>0 /0 ,s 

-0_ n (x) = C_ n , i N e 2 X M(-n, -g + -,x 2 ) ^ . ._„ , . (84) 

' ' T(p + i-n) 2 -\x\~ 9 for x < y J 

with 

e_ „ = 2n - 2g + 1 (85) 

This completes the determination of the spectrum of the supersymmetric pair of Hamiltonians, H + 
and B.-. We note that we have many more states than would be needed from the point of view 
of supersymmetry and, therefore, it is crucial to understand the solutions in a systematic manner 
which we do next. 
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Analysis of the Result 



The large number of solutions obtained is really very interesting and to appreciate their presence, 
let us analyze their behavior in some detail. First, let us note that (as we had pointed out in eq. 
( |65| ) and in the subsequent discussion) a super symmetric ground state would exist for our system 
only if g > — i. Indeed, there exists such a branch of the solutions in our model, namely, the even 
solutions of H+ and one of the two sets of odd solutions of They are, in fact, degenerate in 
energy except for the ground state which is an even state with vanishing energy belonging to the 
spectrum of H + . The solutions have the relative odd parity appropriate to be superpartner states 
and, in fact, it is not hard to check using the properties of the confluent hyper geometric functions 
[17], that, for g > — ^, 



= -v^-^VCn+i^ (86) 

as we would expect from states which are superpartners of each other. This shows that the system 
of solutions, at least, contains the supersymmetric set that we are interested in. 



This, therefore, raises the question about the roles of the other solutions that we have found. 
Both H + and H- have solutions on another common branch, namely, for g < \ and it is worth 
investigating whether, they, too, define a supersymmetric set of solutions. They are degenerate in 



energy (see eqs. (76) and (pT|) ) and they have the correct relative odd parity structure (solutions 



of H + are odd while those for H_ are even). However, it is easily verified that 

Q\g)^{x)^^ n+1 {x) (87) 
even up to normalization constants. On the other hand, it can also be checked equally easily that 

= -v^^^T +1 (*) (88) 
u — ,n+l 

Thus, there appears to be some sort of supersymmetry for these states, but there is no apparent 
state of vanishing energy. Furthermore, there is, of course, the third set of solutions of H + and H_ 
which do not even share the same branch. 



The meaning of all these solutions becomes quite clear once we realize that the supersymmetriza- 
tion of a given Hamiltonian, in this case, is not unique. There are, in fact, various possibilities. For 
example, 

(i) we can choose as the superpotential 

W(x) = --x (89) 
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as we have done so that the supersymmetric pair potentials would be 

V+(x) 
V-{x) 



9 -^ + x 2 -2g-l 
9 -^ + x 2 -2g + l 



(90) 



Here, as we have already seen, a normalizable ground state would exist only for g > — \ (The 
ground state belongs to H + ). 



(ii) If, on the other hand, we choose 



W(x) 



9-1 



— x 



X 



the supersymmetric pair of potentials become 



V+(x) = \ 



V-{x) 



(9-D(9-2) +x2 _ 2g+1 



9 ^ 9 9 - + x 2 - 2g + 3 

x z 



where a normalizable ground state would exist only if g > \. 



(in) Similarly, we could have chosen 



W(x) = - 9 -^-x 



X 



which would lead to a pair of supersymmetric potentials 

V+(x) 
V-{x) 



(fl+l)(0 + 2) n 

Jyy -^+x 2 + 2g + l 



yKh 2 ' +x 2 + 2g + 3 



x^ 



where a normalizable ground state would exist only for g < — \. 



(91) 



(92) 



(93) 



(94) 



In some sense, this is the generalization of the idea that a spin ^ particle, say for example, 
can belong to distinct multiplets of the form (0, ^) or (^, 1). (As is clear now, the transformation 
g <-> — (g + 1) or g <-> —(g — 1) alluded to earlier would simply take one supersymmetric system to 
another and cannot be a symmetry of one supersymmetric system.) 



(iv) However, there is also the possibility of choosing 



W(x) = -i--x 
x 



(95) 
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giving the supersymmetric potentials 



V + (x) 
V-{x) 



9(9 + 1) n 

yyy Z ' + x 2 + 2g-l 



g{g 



+ x 2 + 2g + 1 



(96) 



Namely, H + and H_ could also have reversed their roles (This basically determines which potential 
contains the ground state). In this case, a normalizable ground state would exist only for g < |. 



The meaning of all the solutions is clear now. First, the solutions corresponding to the branch 
g > — 2 define the supersymmetric solutions corresponding to the supersymmetrization in eq. (j§9|). 
The solutions corresponding to the branch g < | also define a supersymmetric set of solutions 
corresponding to the supersymmetrization in eq. (p5|). (It is clear now why the solutions were 
related to each other as in eq. (p8|)-) The fact that the spectrum contains a zero energy state, as 
we would expect in a supersymmetric system, follows from the fact that the potentials in (|9^) are 
shifted with respect to those in (p0[). Finally, the other solutions that do not even share the same 
branch correspond to the supersymmetrizations in eqs. (|9l|) and (|93|). Supersymmetry is manifest 
in spite of the singular behavior of the potentials. It is also clear from this analysis that distinct 
solutions really correspond to distinct supersymmetrizations and in deriving conclusions regarding 
supersymmetry, one should be very careful in identifying the correct solutions from among the 
whole set. 



5 Algebraic Solution 



The supersymmetric system studied in the previous section is a very special system which can also 
be solved algebraically. Therefore, to further understand the properties of this system, we present 
here a short derivation of the algebraic solution of this system. 



Let us consider the supersymmetric system of potentials in eq. (|64|). Defining the supercharges 
as in @, namely, 

1 M s 



we can write the supersymmetric pair of Hamiltonians also as 



H+(g) = QHg)Q(g) 

H-{g) = Q(g)Q\g) (98) 
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However, from the structure of the potentials in eq. (p4|), it is clear that we can write 

H-{g) = Q(g)QHa) = Q ] {g + l)Q(jg + 1) + 2 = H + ( g + l) + 2 (99) 

Such a Hamiltonian, H + (g), (equivalently, the potential) is known as shape invariant [18,12]. The 
ground state of H + {g) satisfies 

Q(g)^P+,o(g, x) = -|= + --x) i)+,o(x) = o 



^J2 \dx x 



or, ip +i0 (g,x) » x 9 e i x (100) 
This solution is normalizable and physical only if 

9>~\ (101) 

Furthermore, the ground state energy of H + (g), in such a case, is zero, 

e + ,o = (102) 

Since H + (g) and H_(g) are supersymmetric partners and share all the energy levels except for 
the ground state of H + (g), it is clear that the ground state energy of H_(g) would correspond to 

the first excited state energy of H + {g) and from the relations in eqs. ( |99|) and ( |102| ), it follows, 
therefore, that the energy of the first excited state of H + (g) is 

e+,i = 2 (103) 



Let us note here that the ground state wavefunction of H_(g) (from the structure in eq. (|9£ 
would follow to have the form 

i; +fi (g + l,x)^x 9+1 e- 1 2 x2 (104) 

and would be normalizable only if 

3 

9> -2 



and this would automatically hold if eq. ( |101| ) is true. 

We can, in fact, construct a sequence of Hamiltonians, in such a case, of the form 

ff<°> = H + ( g ) = QHg)Q(g) 

F« = H^g)=Q(g)Q\g) = Q\g + l)Q(g + l) + 2 
H {2) = Q(g + l)QHg + l) + 2 = Q\g + 2)Q(g + 2) + 4 

H {s) = Q\g + s)Q(g + s) + 2s (105) 
From the structure of the Hamiltonians, it is clear that the ground state of any Hamiltonian 



in the sequence would be normalizable if eq. (101) holds. Furthermore, every adjacent pair of 



Hamiltonians, and H ( s ^ would define a supersymmetric pair of Hamiltonians and share all 
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the energy levels except for the ground state of ij( s— x ) which would have an energy eigenvalue 
2(s — 1). The energy levels of our original pair of Hamiltonians, H + {g) and H_(g) now follow from 
this to be 

e+, n = 2n, e_, n = 2(n + 1) (106) 
These are of course, the energy levels obtained in our earlier analysis (see eqs. (fn]) and fl83[)). 

The wave functions for the original Hamiltonians can also be obtained from an analysis of the 
spectrum of this sequence of Hamiltonians. It is easily seen that [19] 

i>+,n{9, x ) « Q t (#)'0+,n-l(£ + (107) 

which can be iterated to give 

V>+,nG?, x) « QHg)QHg + i)---Q ] {g + n- l)ip+, (g + n, x) (108) 

And since, we know the ground state wave function, V+,o ( see eq. ( |100[ ) ) , all the wave functions 
for the higher states can be explicitly constructed. Furthermore, once we know ip+ >n (g,x), the 
superpartner states can also be obtained from the relation that (see (g) and the discussion following) 

tp-,n(g, x) oc Q{g)i) + , n {g, x) (109) 



It is worth asking if these algebraically constructed states also coincide with the states we have 
constructed in eqs. (72) and (S3) respectively. It is easy to verify the relation that 



d g 

dx x 



x e 



It 2 

z x 



9+1 M(-n + l,g + -,x 2 



(2g + 1) e"^ 2 x 9 M(-n, g + 



1 



(110) 



which is the appropriate recursion relation in eq. (|1Q7| ) for our system. Furthermore, noting the 



form of the ground state wave function in eq. (100) 



x 9 



1„2 
2 X 



x 9 M(0,g + -,x 2 



(111) 



and choosing an even function for the ground state, as we normally do, (as far as we can see, the 
algebraic method cannot determine a priori whether the function should be even/odd), it is clear 
that all the states associated with H + (g) would be even states and would coincide with those in 
eq. (|72|). Once these are identified, the superpartner states also follow to coincide with those in 
eq. (^) and are odd. This concludes the algebraic construction of the supersymmetric spectrum 
of the system and shows again that supersymmetry is manifest. 



6 The Puzzle 



An interesting puzzle has been raised in the literature [11] in connection with the supersymmetric 
system of eq. ( |64[ ) and in this section, we describe the resolution of the puzzle from a systematic 
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analysis of the problem. Namely, let us look at the pair of supersymmetric potentials in eq. (64) 



V+(x) 
V-{x) 



9(9 ~ 1) 



2g-l 
2g + l 



and note that, for g = 1, the potentials take the form 

V+(x) = 

V-¥> = + 



1 



(112) 



As we have noted earlier, g = 1 is not particularly a special point in the parameter space (only 
g = is) and, consequently, our analysis would continue to hold for this value of the coupling as 
well. Furthermore, this value of the coupling is consistent with g > — ^ and, therefore, we would 
expect the spectrum to be supersymmetric with the ground state energy of H + (g = 1) vanishing as 
the analysis of the last two sections have shown. On the other hand, if we look at the potential V+ 
in eq. (|l!2j ), it is clear that it is a harmonic oscillator potential with the zero point energy shifted 
by — |. It would appear, therefore, that the ground state energy of H + (g = 1) would be negative 
contrary to what our analysis has shown and in violation of the theorems of supersymmetry (The 
supersymmetry algebra tells us that the ground state energy in a supersymmetric theory has to 
be positive semi-definite.). This, therefore, raises an interesting puzzle and it has been discussed 
in the literature from various points of view. However, the resolution of this puzzle is really quite 
simple as we show next. 



In concluding that the potential V+ in eq. (112) would lead to a negative energy, we have, of 



course, tacitly assumed the ground state wave function for the system to be the standard oscillator 
wavefunction, namely, 

i/>o(x) oc (113) 

which, as it should be, is an even function. On the other hand, from our analysis of the last two 
sections (see eq. (|72|) or (|100|) ), we would conclude that the ground state wave function for our 
supersymmetric system would have the form 

1 2 

if> +7 o(x) oc e~z x \x\ (114) 

It can also be explicitly checked that this indeed has a vanishing energy consistent with our general 
analysis. Therefore, understanding the question of the negative energy reduces to understanding the 
role of the wave function in eq. ( |113| ). It is, of course, clear from the explicit spectrum of the system 
we have constructed that it is not part of our supersymmetric system. It must, therefore, belong 
to a different supersymmetric system since as we have seen earlier, distinct supersymmetrizations 
lead to distinct forms of the wave functions. 



In fact, let us choose 



W{x) = -- — -- x (115) 
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which would lead to the pair of super symmetric potentials of the forms 



V+(x) 
V-(x) 



9(9 ~ 1) 



+ x 2 + 2g - 3 



(116) 



It is clear that this supersymmetric system would have a normalizable, even ground state of the 
form 



ipo(x) oc \x 



1-9 



e 2" 



(117) 



which is normalizable for g < |. It is also clear that this wave function coincides with ( |113|) when 
5 = 1 (which is in its range of validity). In other words, the wave function in eq. (|113[) does 
not belong to the supersymmetric spectrum of (|64|), rather, it is a legitimate ground state of the 
supersymmetrization in eq. ( 115 ). For g = 1 the potentials of (116), of course, have the form 



V+(x) 
VL(x) 



1 r 

2 
1 

2 



x 



x 2 + 1 



(118) 



which is nothing other than the supersymmetric oscillator which, as we know, has a zero energy 
ground state. 



Thus, there is really no puzzle and the theorems of supersymmetry are, indeed, on firm footing. 
The crucial thing to learn from this analysis is what we have emphasized earlier, namely, in dealing 
with supersymmetry in such a system it is essential to identify the appropriate solution relevant 
for a particular supersymmetric system because every distinct solution corresponds to a distinct 
supersymmetrization of the system. 



7 Conclusion 



In this paper, we have studied systematically two classes of supersymmetric quantum mechanical 
models - one consisting of a singular boundary and the other with a singular potential. We have 
shown that, contrary to the conventional understanding [10-12], supersymmetry is manifest in these 
systems. In particular, for a system with a singular potential such as -7, we have shown that the 
solution of the Schrodinger equation leads to several distinct solutions corresponding to distinct 
supersymmetrizations of the system. Consequently, it becomes quite important to identify the 
appropriate wavefunctions when supersymmetric properties are being investigated. We have also 
solved the system with a singular potential algebraically using the ideas of shape invariance which 
further corroborates the results of our direct analysis of the Schrodinger equation, namely, that 
supersymmetry is manifest in this system. Furthermore, we have shown how a puzzle raised in 
the literature finds a natural resolution when one does a careful analysis restricting oneself to the 
appropriate branch of the supersymmetric solutions. Finally, we would like to conclude by noting 
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that supersymmetry is known to be robust at short distances (high energies). The singularities 
discussed in the quantum mechanical models occur at short distances and, therefore, it is intuitively 
quite clear that they are unlikely to break supersymmetry. Our detailed, systematic analysis only 
reinforces this. 
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